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We report on measurements of differential cross sections da/dpx for prompt charm 
meson production in pp collisions at y/s = 1.96 TeV using 5.8±0.3pb _1 of data from 
the CDF II detector at the Fermilab Tevatron. The data are collected with a new 
trigger that is sensitive to the long lifetime of hadrons containing heavy flavor. The 
charm meson cross sections are measured in the central rapidity region |y| < 1 in four 
fully reconstructed decay modes: D° — ► K~ir + , D* + — > D°ir + , D + — > K~tt + it + , 
Df — ► 07r + , and their charge conjugates. The measured cross sections are compared 
to theoretical calculations. 

PACS numbers: 12.38.Qk,13.85.Ni,13.25.Ft,14.40.Lb 

Measurements of the production cross sections of hadrons containing b quarks or 
charm quarks (heavy flavor hadrons) in pp collisions provide an opportunity to test pre- 
dictions based on Quantum Chromodynamics (QCD). Previous measurements of B meson 
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production cross sections in pp collisions at y/s = 1.8 TeV were about three times larger 
than Next-to-Leading-Order (NLO) QCD predictions Q], although recent calculations with 
a more accurate description of b quark fragmentation have reduced this discrepancy to a 
factor 1.7 jj, 0]. Charm meson production cross sections have not been measured in pp 
collisions and may help with understanding this disagreement. The upgraded Collider De- 
tector at Fermilab (CDF II) has a new capability to trigger on tracks displaced from the 
beamline originating from the decay of long-lived hadrons containing heavy flavor quarks. 
We report here measurements of prompt charm meson cross sections using data recorded 
with this trigger in February and March 2002, corresponding to 5.8 ± 0.3 pb -1 of integrated 
luminosity. 

r"i 

An overview of the CDF II detector can be found elsewhere p], only the components rele- 
vant to this analysis are described here. The CDF coordinate system has the z axis pointing 
along the proton momentum; ip is the azimuth, 9 is the polar angle, and r is the distance 
from the proton beam axis. The CDF II central tracking region covers the pseudorapidity 
region \t}\ < 1, where rj = — ln[tan(#/2)]. A superconducting solenoid provides a nearly 
uniform axial field of 1.4 T. The silicon vertex detector (SVXII) |7j consists of double-sided 
microstrip sensors arranged in five concentric cylindrical shells with radii between 2.5 and 
10.6cm. Surrounding the SVXII is the Central Outer Tracker (COT) ^|, an open cell 
drift chamber covering radii from 40 to 137 cm. The COT has 96 layers, organized in 8 
superlayers, alternating between axial and ±2° stereo readout. 

CDF II has a three-level trigger system. We describe here the trigger used in this analysis. 
At the first trigger level, charged tracks are reconstructed in the COT axial projection by a 
hardware processor (XFT) 9]. The trigger for hadronic charm decays requires two oppositely 
charged tracks with px > 2GeV/c and the scalar sum of the p^s larger than 5.5GeV/c, 
where pt is the magnitude of the component of the momentum transverse to the beam 
axis. At the second trigger level, the Silicon Vertex Tracker (SVT) [l^ associates axial strip 
clusters from the four inner SVX II layers with XFT track information. The SVT measures 
the distance of closest approach of a track relative to the beam axis (impact parameter or 
do) with a resolution of 50 /mi, which includes a contribution of 30 fim from the beam spot 
transverse size. Events containing hadronic decays of heavy flavor hadrons are selected by 
requiring two tracks with 120 fim < d < 1 mm each. At the third trigger level, a farm of 
computers performs complete event reconstruction online; the opening angle Sip between the 
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two trigger tracks is required to be between 2° and 90°, and the intersection point in the 
r-(p plane projected along the net momentum vector of the two tracks must be more than 
200 ^m from the beamline. 

We reconstruct charm mesons in the following decay modes: D° —>■ K~tt + , D* + — ► D tt + 
with D° — > K~tt + , D + — > K~tt + 7t + , — > 07r + with — > K + K~ , and their charge conju- 
gates. For every track pair that satisfies the trigger requirements (trigger pair), we form 
one D° — > i^~7r + candidate and a second candidate with the mass assignments swapped. 
No particle identification is used in this analysis. D° candidates within three standard de- 
viations of the D° mass are combined with a third track with pt > 0.5GeV/c to form 
D* + — > D°7T + candidates. The three-body decays of the D + and are reconstructed by 
combining a trigger pair with a third track having axial hits in at least three out of five 
SVXII layers and performing a vertex fit based on axial track information only. For Df 
reconstruction, we specifically require the K~ir + pair to satisfy the trigger requirements, 
since the typical opening angle between two kaons from <p decay is close to the 5ip > 2° 
trigger requirement. Each <p candidate is required to have a mass within ±20MeV of the 
world average mass . The D meson candidates are binned in px as indicated in Table HI 
The signals summed over all pt bins are shown in Fig. ^ 

The D° yield is obtained from a binned maximum likelihood fit to the K~ir + invariant 
mass distribution, with a linear function for the combinatoric background, a narrow Gaussian 
for the D° signal, and a wide Gaussian with the same normalization describing D° — > K~tt + 
with the wrong mass assignment. We determine the shape of the mass distribution resulting 
from the wrong mass assignment using D° from D* + decay where the charge of the low- 
momentum pion determines the mass assignment of the D° decay products. The D* + yield 
is extracted from the distribution of Am = m(K~ix + 'K + ) — m(K~7T + ), the mass difference 
between the K~tc + it + and the K~7t + combination. The signal is modeled with two Gaussians 



with equal means, and the background is characterized as a a/ Am — m„-exp(6(Am — m,)), 
where a and b are free parameters in the fit. The D + signal is described with two Gaussians 
and the background is described with a linear function. In the 07r + mode, we model the 
invariant mass distribution as a linear background with two Gaussians, one corresponding 
to the D + and one to the D+. We find 36804 ± 409 D° -^K~7r + , 5515 ± 85 D*+ -> D°n + , 
28361 ± 294 D + -> K-n + n + , and 851 ± 43 D+ -> 0tt + , where the uncertainties quoted 
are statistical only. We vary the signal and background models and attribute systematic 
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uncertainties on the signal yield in the range of 1% — 6%, depending on the decay mode and 
the pt range of the candidates. 

We separate charm directly produced in pp interactions (prompt charm) from charm 
from B decay (secondary charm) using the impact parameter of the net momentum vector 



of the charm candidate to the beamline [12J. Prompt charm mesons point to the beamline. 
The shape of the impact parameter distribution of secondary charm is obtained from a 
generator- level NLO Monte Carlo (MC) simulation of B meson production Q and decay , 
smeared with a resolution function (Gaussian + exponential tails) obtained from a sample of 
Kg — > 7r + 7r~ decays that satisfy the trigger requirements. The impact parameter distribution 
of the reconstructed charm samples, shown for the D° in Fig. |2J is fit to a prompt and a 
secondary component. The prompt fraction is measured for each p? bin. Averaged over all 
p T bins, (86.6 ± 0.4)% of the D° mesons, (88.1 ± 1.1)% of D*+, (89.1 ± 0.4)% of D+, and 
(77.3 ± 3.8)% of are promptly produced (statistical uncertainties only). The systematic 
uncertainties on the prompt fractions are estimated by removing the non-Gaussian tail in 
the resolution function and evaluating the variation. The relative uncertainty is found to be 
in the 3% — 4% range, depending on the decay mode. 

Using a hit-level simulation of the COT, overlaid with data events from the hadronic heavy 
flavor trigger to reproduce a realistic occupancy, we measure a reconstruction efficiency in 
the COT of 99% for tracks with p T > 2 GeV/c, falling to 95% at p T = 0.5 GeV/c. The 
efficiency for finding three SVX II axial hits on a reconstructed track is measured from data 
to be about 85%. The efficiencies of the trigger hardware XFT and SVT to reconstruct 
tracks are measured from data samples without XFT or SVT requirements [l2|. The XFT 
tracking efficiency is greater than 95%. In the data-taking period considered, the SVX II and 
the SVT were not yet fully operational, and the efficiency varied as certain SVX II modules 
were included or excluded from the trigger. Therefore, we measured the SVT efficiency in 42 
periods, each corresponding to one pp store of the Tevatron, and characterized the efficiency 
as a function of the track azimuth (p, the longitudinal position zq, the polar angle 9 and the 
transverse momentum p?. The average single-track efficiency of the SVT for this early data 
taking period was about 42%. 

The measured efficiencies are applied to a generator-level NLO MC simulation of charm 
meson production and decay to calculate the trigger and reconstruction efficiencies, taking 
into account decay in flight and hadronic interactions of the charm meson decay products. 
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The MC pr spectrum of the charm mesons is reweighted to match the measured px spectrum. 
The integrated cross section cij in each px bin i with \y\ < 1 (where y = | ln( ^^ z ) and E 
is the energy of the charm meson) is calculated using the following equation: 

N*/2 ■ f D>i 



(J; 



(1) 



f£dt-ei-B' 

where iVj is the number of charm mesons in each px bin and fo,i is the fraction of prompt 
charm in that bin. The integrated luminosity / Cdt at CDF is normalized to an inelastic cross 
section of a vS = 60.7±2.4mb llJl. The rate of inelastic colisions is measured with Cherenkov 
luminosity counters [15[ and has an uncertainty of 4.4%. The factor | is included because we 
count both D and D mesons, but we report cross sections for D alone. We verified that the 
D and D cross sections are equal within statistical uncertainties. The branching fractions 

n 

B are taken from Ref. pjj. For the D cross section, we sum the branching fractions of 
D° — > K~n + and D° — > K + tc~ , since both contribute to the observed signal. The combined 
reconstruction and trigger efficiency q varies from 0.12% to 1.9% depending on the decay 
mode and the px bin. Systematic uncertainties on the trigger and reconstruction efficiency 
arise predominantly from the uncertainty on single-track efficiencies and two-track efficiency 
correlations. They also have contributions from ionization energy loss, hadronic interactions 
in the inner tracker material and the size of the interaction region. The combined systematic 
uncertainty on the trigger and reconstruction efficiency is in the range of 8%— 14%, depending 
on the decay mode and the px range of the D mesons. 

The total cross sections are obtained by summing over all px bins. However, the last px 
bin is replaced by an inclusive bin with px > 12 GeV/c. We find a(D°,px > 5.5 GeV/c, \y\ < 
1) = 13.3 ± 0.2 ± 1.5 fib, a(D* + ,p T > 6.0GeV/c, \y\ < 1) = 5.2 ± 0.1 ± 0.8 /A, a(D + ,p T > 
6.0GeV/c, \y\ < 1) = 4.3 ± 0.1 ± 0.7 //b and a{D+,p T > 8.0GeV/c, \y\ < 1) = 0.75 ± 0.05 ± 
0.22 fib, where the first uncertainty is statistical and the second systematic. To calculate the 
differential cross sections, we divide Oi by the width of the px bin. Since we report da /dp T 
at the center of each px bin, we apply a correction to account for the non-linear shape of 
the cross section, using the px reweighted MC to obtain the shape of the cross-section inside 
each px bin. The results are listed in Table HI 

The measured differential cross sections are compared to two recent calculations , 
as shown in Fig. El The uncertainties on the calculated cross sections are evaluated by 
varying independently the renormalization and factorization scales between 0.5 and 2 times 



10 



the default scale. Ref. [16 uses a default scale of ym^ + p?, where m c = 1.5GeV/c 2 is the 
c quark mass, while ReTQ uses a default scale of 2 J ml + p\. Contributions from other 



sources, such as the charm quark mass, the value of the strong coupling constant and the 
fragmentation functions, were reported to be smaller and are not taken into account. 

In conclusion, the measured differential cross sections are higher than the theoretical 
predictions by about 100% at low p? and 50% at high p?. However, they are compatible 
within uncertainties. The same models also underestimate B meson production at ^/s = 
1.8 TeV by similar factors 3,0,3- 
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da(\v\ < D/dvr 


[nb/(GeV/c)l 






(jfiirfcrs! T) r r 


D° 


D* + 


D+ 


D+ 


fGeV/cl 


[GeV /cl 










5.5 — 6 


5.75 


7837 ± 220 ± 884 








6-7 


6.5 


4056 ± 93 ± 441 


2421 ± 108 ± 424 


1961 ± 69 ± 332 




7-8 


7.5 


2052 ± 58 ± 227 


1147 ± 48 ± 145 


986 ± 28 ± 156 




8-10 


9.0 


890 ± 25 ± 107 


427 ± 16 ± 54 


375 ± 9 ± 62 


236 ± 20 ± 67 


10 - 12 


11.0 


327 ± 15 ±41 


148 ± 8 ± 18 


136 ± 4 ± 24 


64 ± 9 ± 19 


12 - 20 


16.0 


39.9 ±2.3 ±5.3 


23.8 ±1.3 ±3.2 


19.0 ±0.6 ±3.2 


9.0 ± 1.2 ±2.7 



TABLE I: Summary of the measured prompt charm meson differential cross sections and their 
uncertainties at the center of each px bin. The first error is statistical and the second systematic. 
The products of the branching fractions Q used are (3.81 ± 0.09)%, (2.57 ± 0.06)%, (9.1 ± 0.6)% 
and (1.8 ± 0.5)% for LP, D*+, D + and D+, respectively. 




pi I pi i i i I 

1.7 1.8 1.9 2 1.8 1.9 2 2.1 

m(K "it) [GeV/c 2 ] m(K*K it*) [GeV/c 2 ] 



FIG. 1: Charm signals summed over all px bins: (a) invariant mass distribution of D — > K~ir + 
candidates; (b) mass difference distribution of D* + — > L>°7r + candidates; (c) invariant mass dis- 
tribution of D + — ► K~ir + ir + candidates; (d) invariant mass distribution of D + — > (J)tt + and 
Df — > <j)7r + candidates. The curves show the results of the fits described in the text. 
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Impact Parameter [|im] 

FIG. 2: The impact parameter distributions of the D° mesons, measured from the ±2<r signal 
region of the invariant mass distribution and corrected for combinatoric background measured in 
the invariant mass sidebands. The solid curve is the fit result summed over all pt bins. The dashed 
curve shows the contribution of secondary charm from B decay. 
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FIG. 3: The measured differential cross section measurements for \y\ < 1, shown by the points. 
The inner bars represent the statistical uncertainties; the outer bars are the quadratic sums of 
the statistical and systematic uncertainties. The solid curves are the theoretical predictions from 
Cacciari and Nason 16|, with the uncertainties indicated by the shaded bands. The dashed curve 
shown with the D* + cross section is the theoretical prediction from Kniehl rj|]; the dotted lines 
indicate the uncertainty. No prediction is available yet for production. 
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